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ABSTRACT 

We study the intermittency properties of two branching processes, 
one with a uniform and another with a singular splitting kernel. The 
asymptotic intermittency indices, as well as the leading corrections to 
the asymptotic linear regime are explicitly computed in an analytic 
framework. Both models are found to possess a monofractal spectrum 
with ip q = q — 1. Relations with previous results are discussed. 
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1 Introduction. 



The original proposal of intermittency in multiparticle production con- 
sisted of a method of quantifying and measuring short-range correlations 
in multi-hadron final states, as well as a hypothesis about the self-similar 
nature of fluctuations. This fractal character is supposed to have dynamical 
origins and expected to provide valuable information about the underlying 
mechanism of hadron production. The subject of intermittency has been 
intensely pursued since its introduction in Q, both theoretically [§, ||| and 
experimentally [||, ||] and is by now a standard topic in soft hadronic physics. 

On the theoretical side, the property of intermittency was found to hold 
in many models, most remarkably a-models Q, of which detailed studies 
have been made ([§,0] an d references therein). A different class of systems 
where intermittency is expected to be present is branching processes. These 
have been applied to the phenomenological description of high-energy pro- 
cesses both of electromagnetic Q and hadronic || 10, |IT) nature. In the 



context of intermittency, branching processes provide a laboratory to study 
short-range correlations and fluctuations in multi-particle production either 
analytically or numerically. 

Intermittency in branching processes has been studied in [[12], where 
a proof of intermittent behavior is offered. It has also been considered 
in the more specific framework of QCD branching in [13, 14 1, where the 
problem of infrared divergences arises. In this paper we adopt a point of 



view complementary, in some sense, to that of |13|, 14] and closer to [12]. 
We shall study the intermittent regime of two one-species branching models, 
characterized by a uniform and a singular splitting kernel, respectively. Our 
aim is to determine whether these mathematical models display intermittent 
behavior, of what kind, in what limit, and with what type of corrections. 
We are also interested in the relationship [15] between intermittency and 
KNO scaling 16]. Our mathematical approach is different and, we believe, 
simpler than that of [12], which allows us to obtain more specific results. 
The price to pay for those results is that we have to give up generality, by 
considering two particular cases, and that we defer for future consideration 
the problem of infrared-divergent kernels. 

The outline of the paper is as follows. In Section 2, we define the model, 
explain those general features of branching processes which are relevant to 
our purposes, and fix our notations and conventions. In Section 3, we study 
the particular models in terms of evolution equations for inclusive distri- 
butions. For both splitting kernels, a monofractal spectrum with maximal 
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intermittency indices [17] is found, as well as the leading corrections to the 
asymptotic behavior of scaled factorial moments and a precise characteriza- 
tion of the regime where intermittency appears. In section 4 , a summary 
of results is given, together with our final remarks. An Appendix, finally, 
gathers some complementary material related to the mathematical approach 
used in the main body of the paper. 



2 Model. Notations and Conventions. 

The branching process under consideration consists of particles characterized 
by an energy fraction x and a virtuality t, the latter being the evolution 
parameter of the system. Mathematically it is defined by the splitting kernel 
P(x) entering the evolution equations for the probability V to having a 
branching at virtuality t, 



dP 
dt 



dxP(x) (1) 



This equation is inspired in the Altarelli-Parisi-Gribov-Lipatov equations 
for partonic evolution |l8|, [l9|], with the simplificatory assumption that the 
running coupling constant a s (Q 2 ) does not depend on x, Q 2 = Q 2 (t), so 
that it can be absorbed in the evolution variable t In this way, a 

stationary Markov branching process is defined and evolution equations for 
the generating functional of transition probabilities can be found. We shall 
omit the derivation here, which is given in detail in, e.g., |2^, [2l], ^] (see 
also the Appendix, where the relevant equations are summarized). 

We shall be interested, as will be made clear later in this section, in 
inclusive transition probabilities D n (x\, . . . ,x n ;t) which represent the prob- 
ability to observing particles with energy fractions x%, . . . ,x n , as well as 
particles with any other values of x, at virtuality t, and starting with one 
particle with x = 1 at t = 0. These inclusive distributions are normalized 
to the factorial moments of multiplicity, 

/ dxi . . . dxkD^xx, ...,x k ;t) = (n(n - 1) • • • (n - k + (2) 
Jo 

where the multiplicity n is the total number of particles in the system at vir- 
tuality t. This normalization condition follows naturally from the definition 
of D n in terms of the generating functional of the process (see Appendix 
and Igl,!!). 
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We shall further assume that the spliting kernel P(x) is normalized to 
unity, which implies no loss of generality in the one-species case; that it is 
symmetric (P(x) = P(l — x)), which follows from energy conservation; and 
that its support is contained in the interval [0, 1]. Moreover, we shall adopt 
the convention that all probability densities, inclusive or exclusive of any 
order, vanish outside [0, 1]. 

The inclusive densities of order n, D n (xi, . . . ,x n ,t), satisfy the forward 
equation (see Appendix), 

8D n 



dt 



- (x\ , . . . , X n , t) 



-nD n (xi, ... ,x n ,t) + 2 I dz P(z) ^ -D n (xi, x n , t) 

Jo k=1 z z 



k>j=i 



+ 2 p [ — T — — T — D n-l(xi,...,Xj+X k ,...,Xk,...,X n ,t)(3) 

Xj ~r %h I Xj ~~\- Xfc v v y 
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where the hat indicates that the variable Xk is omitted and we used the 
symmetry of P{z). The initial condition is given by, 

D n (xi, . . . , x n ,t = 0) = 5 ln 5{x 1 - 1) (4) 

We now define the factorial moment densities of order n, (cf., eq. (|2|)) 

d n (x,t) = [ dx 2 - ■ ■ dx n D n (x,x 2 , ■ ■ ■ ,x n ,t) (5) 
jo 

Separating x\ in the equation for D n and integrating dx 2 ■ ■ ■ dx n we obtain 
the forward evolution equation for d n , 

dd 

-Qf(x,t) = (n-2)d n (x,t) + (n-l)(n-2)d n „ 1 (x,t) 

+ 2 [dz P(z)-d n (-) + 2(n - 1) [dz P(z)~d n -i (-) (6) 
Jo z \zj Jo z \zj 

with initial condition, 

d n (x,t = 0) = 5 ln 5(x-l) (7) 
The Mellin transformed equation, 

" a 's,t) = (n-2 + 2P(s))(d n {s,t) + (n-l)d n - 1 (s,t)) (8) 



dt 

d n (s,t = 0) = 5 ln (9) 
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can be explicitly solved, to obtain, 

, , / vn-i ~, T(n-l + 2P(s)) 

d n ( s ,t) = e ("- 2 )* fl - e"*) n e 2P W ^ . , J (10) 
V ' r (2P( 8 )) 

Substituting s = 1 in this expression, we obtain the factorial moments cor- 
responding to a pure-birth, binary-fission process, which are independent of 
the form of the splitting kernel p^j, 

(n(n - 1) • • • (n - q + 1)) = g!e 9i (l - e"*) 9 " 1 (11) 

We shall find useful later the following expanded form of the solution, 

_1 _ n— 1 

£(.,t) = (l - e-*)"" 1 e 2P ^ £ C^" 1 (2P(s))' 

where 



121 



(13) 



[Cm = Z)^>...>i r =l *1 " " " * 

with r = n — 1 — m and n > 2 



2.1 Intermittency indices. 

Our goal in this paper is to obtain the intermittency indices of a branching 
process, defined as the logarithmic slope of reduced factorial moments as 
functions of bin size Q. Reduced factorial moments are defined as, 

F.= | (14) 

U q = (n ■ ■ ■ (n — q + 1)) being the factorial moment of the multiplicity within 
the given bin. Notice that factorial moments U q are obtained from the above- 
defined densities as, 

U q (x;t) = J dx dg(x;t) (15) 

where the integration extends to the bin under consideration. 

Due to the asymmetry of the process, (x)(t) being a monotone decreasing 
function of t (see Appendix), the relevant bin in any partition of [0,1] is 
the one adjacent to the point x = 0. For that reason, it is appropriate to 
apply vertical analysis (see, e.g., [|], |24|) to that bin in order to find the 
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intermittency slopes. However, the particle density U\ need not be uniform 
within the bin under consideration, so that we shall define intermittency 
indices in terms of F q ELS £1 function of U\ [Eq], 

ate 

This is the definition of intermittency indices that we shall adopt in the 
sequel. 

Before going to a particular case in the next section, we would like to 
mention another consequence of the asymetry of the process. For long enough 
t, and a fixed resolution in x, almost all particles will be contained in the 
bin adjacent to x = 0, where the process will continue its evolution as a 
usual branching process with only multiplicity degrees of freedom. But this 



kind of processes KNO scale [10, 16, leading to constant reduced factorial 
moments in the limit t — > oo, and consequently to zero intermittency indices 
in that limit. We shall explicity verify this fact below . 



3 Two Particular Cases. 

In this section, we shall compute the asymptotic form of reduced factorial 
moments with the expressions derived in the previous section. Factorial mo- 
ment densities can be exactly found for P(x) = 1 and P{x) = 8{x — 1/2), 
which we shall exploit to obtain the asymptotic moments at a given resolu- 
tion. 



3.1 Uniform Splitting Kernel. 

For the "uniform model", we have P{z) = 1, P(s) = 1/s. By expanding the 
exponential in eq. (12), we find 

(2t) n 1 



di(s,t) 



E 



n=0 



n! s n 



(17) 
00 (1t\ n 1 

d 2 (s,t) = 2(l-e-*)E^-^ (18) 

71=0 

d n (s,t) = e^Hl-e-r- 1 E Cr^E^TTTO ( 19 ) 



m=l 



fc=0 
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Applying the inverse Mellin transform term by term leads to, 



dx(x,t) = e~\ h (\/-8tlnx) +e~ t 8(x - 1) (20) 

y — m x \ ' 

d 2 (x,t) = 2(1 - e -t )Io (\/-8tlna;) (21) 
d»(z,t) = 2e (n - 2) *(l - e"*)"- 1 x 



n-l 



m—1 



x E^h/^r^l /«-i(v=8*Ei) (22) 



m=l 



where I n denotes a modified Bessel function of the first kind and order n 
|26f| . Since these expressions are given by finite sums, in order to check that 
they are solutions to eq. (||) it is enough to see that their Mellin transform is 
correct. This is easily done by change of variables, turning one-sided-Mellin 
into Laplace transforms and using the known transforms of Bessel functions 

& 

The factorial moments U n (x,t) = Jq dz d n (z,t) are then given by 



U n (x,t) = 2e (n - 2)t (l-e- t ) n - 1 x 

n-l roo / pJT\ m ~ l 

^Tcr 1 / dye-y(J^) I m -l(V&ty) (23) 
m=l J - lnx \ V 1 ) 



X 
m= 



and analogously for U\, U2- We shall now assume that the argument of 
the Bessel functions is large, so that we can approximate them by their 
asymptotic expansion [26] to obtain, 



" 1 poo 

U n (x,t) = 2e("- 2 )*(l - e-'r- 1 ^ C n m x \ dy 

m=l J - lnx 



2y\ m 1 ( 4(m-l) 2 -l 



2vr^/8ty 



1 - — 7= + • • • (24) 



where we have retained the next-to-leading term in order to estimate the 
order of magnitude of corrections later. 



3.1.1 Large t. 

We consider now a resolution x as small as desired but fixed, and t — > oo. 
We can then drop the correction term in eq. (|24|) and re-write the leading 
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term by means of the substitution rj = (y/2t) 1 / 4 as, 



n-1 [7JT 

U n (x,t) = 2e nt (l-e~ t ) n - 1 £ C£-W- 

m=l 



'(tt) 

For i — > oo we can estimate this expression by saddle point in the form, 

n-1 njr 

U n (x,t) 2e n *(l - e"*)™- 1 V C^^J — x 

' V 7T 

m=l 

^2 



x /, ,^e-^il 2m (26) 
or, for p = VSt(i] — 1), 

n-1 iy 

U n (x,t) « 2e"*(l-e-*r- 1 £ C^^Y" x 



x rdpe~ p2 (-P— + \) m (27) 



OO 



2V2t 



Notice that we substituted — oo for — \/8t in the lower limit of the inte- 
gral. Applying the binomial expansion to the bracket in the integrand and 
integrating term by term, the resulting sum shows that in the limit the in- 
tegral takes the value ^/k. Then, by definition of C^ _1 , eq. (13), we find 
U n (x,t) ~ n\e nt , which corresponds to the KNO value given by, 

^ = n\ (t-oo) (28) 

As discussed in the previous section, this leads to zero intermittency indices. 



3.1.2 Small x. 

For t fixed, we now consider the regime of small x. Although the scale of 
smallness will be apparent as we proceed with the computation of factorial 
moments, we can in principle assume x <C (x), where (x) = t/{e t — 1) is the 
average energy fraction (see Appendix). 
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By changing the integration variable in eq. (24) to r\ = ^fy — y2t, we 
obtain, 



n— 1 



U n (x,t) = e^l-e"*) 



-i\n-l 



m — 1 

2\ — 



m-1/2 



x 1 



(m — 1) 



4V2t V (l + ^f) 



+ 



(29) 



where we denoted a = V '— E x — y2t for brevity. By applying the series 
expansion, valid for a > y2t, 



2t\ 



p/2 



E 



1 Tip/2 + 1) 



^ 2n\ r(p/2 - n + 1 
to the previous expression of f7 n , we have 

4 



a] r(?-? + i 1 a 



(30) 



n-l 



I7 n (x,t) = e^l-e" 4 ) 



-t\n-l 



27r(8t)V 



174 E G . 



n-l 



m=l 



TTl — 1 

2\ — 



E 



(V2?) r(m + l/2) 



. , m k 1 



k J 2k\ T(m + 1/2 -k) \ 2 2 4' 
(m-l) 2 -l k 



1 



H 



8t m-1/2 
Some remarks about this expression are in order here, 



(31) 



i. The correction factor in brackets comes from the asymptotic expansion of 
the modified Bessel function. Since it is linear in k, it does not affect the 
rate of convergence of the series as given, e.g., by the ratio test. 

ii. This factor does not affect the zeroth-order term k = 0. Similarly, the next 
order term would not affect the k = and k = 1 terms, etc. 
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iii. As we shall see, only the leading term in the series need be kept. If we 
wanted to retain more terms in the series, we would have to take into ac- 
count the next corrective terms coming from the asymptotic expansion of 
the Bessel function. Such expansion is, however, divergent, so that only a 
finite number of terms may be kept. 

Regarding the convergence of the series, we can use the following asymptotic 
expression, obtained by a saddle-point argument for large uj, 

e~ x x~" 



r(— u, x) 



1 + 



UJ + X + 1 

to obtain the rate of convergence of the series as 



(uj 



+00 



(32) 



(33) 



V-lnx - \[2t 

We see, then, that the series is absolutely convergent for 8t < — In x and, 
furthermore, that for \/2t <C y/— hxx it is justified to truncate it. We shall 
retain the leading (k = 0) and next-to-leading (k = 1) terms, the second one 
in order to have an estimate of the corrections. 



U n (x,t) 

e nt (l - e"*) 

n-l 



n-1 



27r(8i)V4 



X 



Er~m— 1 



m=l 



1 + 



m — 1 

~ 1. 



m 1 



^ r(f-i,a 2 )/ 



(m - if 



1 



1 



8/ 



( 2 1 

We now use the asymptotic form [^] for uj > fixed, x 

1 



m - 1/2 



00, 



T(uj, x) 



1 + 



UJ 



+ 



;s4) 



(35) 



to obtain 

U n (x,t) 



(!-e-) 



n-l 



n-l 

x c*m 1 ( T 



m=l 



m— 1 
— 1 



l + (m-l/2) 



2vr(8t) 1 /4 

1 ^-a^m-3/2^ x 

(m-1) 2 - 



1 



8f 



m - 1/2 



(36) 
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or, recalling that a = \J — In a: — y/2t, 



U n {x,t) ^ e^-V* (l - e"*) 



n-i 2e 



2vr(8t) 1 /4 



1^ ^ m— 1 

m=l 



1 + 



24 



lnx 



(m - l) 2 - 1 
8* 



(37) 



The leading contribution in the limit x — > comes from the m = n — 1 term. 
Keeping only this and the next-to-leading one we finally have, 



U n (x,t) 



>-2)t 



1 + 



1 — e ) /2\ 2 



(8t)V4 



(e^^x(-lnx)f-f) 



24 



lnx 



1 - 



n-2) 2 -l (n-l)(n-2) /2 



84 



+ 



(38) 



Mean value U\ and factorial moment Ui- For the cases n = 1, 2 we 
have, 



EMM) = e~* / dxh(\/-8t In: 
jo 



24 



— lnx 



EMM) = 2(1 -e _t ) /dx/ (V-8tlnx) 
j o 

Expanding the Bessel functions as before, we arrive at 



EMM) 



e*(8i)V4 « 1 (^24) n „/l 



v 7 ^ 



E 



n=0 



2n!r(l/2-n) V4 2 



n 



r ---,a' 



3n 



(l-eV* ^ 1 (v^)" ^3 n 2 
EMM) = ToTvTTi E~ rr(---,a 



(24)V4 ^ Q 2n!r(3/2 



n) 



2' 



1 + 



n 

164 



(39) 
(40) 

(41) 
(42) 



Retaining only the leading and next-to-leading terms and expanding the 
gamma functions results in, 



EMM) 



e~*(84) 1 / 4 e ^- 8tlnx x 
2^ (-lnx) 3 / 4 



1 + 



24 



lnx V 44 



(43) 
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I 1 

X I — I 

\ 2 16t 



(44) 



Intermittency indices. Applying the expression for the intermittency 
indices to the asymptotic expansions found above, we have, 

^ = 1 + ^-FlnW(^ i ^ + I^) + --- ^ 

(n - 1) A , . 

<^ n = (n - 1 + 9 , 1 \ ~ r+l : > 3/2 + • • • 46 
2(-lnaj) Vt(-lnx) 3 / 2 

n+l , — 

A = (n - l)(n - 2 t 2 ~ n/2 - ± — - — (n 2 - 4n - 3) 

We see, then, that for \/2t <C \/ — hix and n <C — In x we have intermittency 
with a monofractal spectrum of indices, with logarithmic corrections that, 
to first order, do not depend on t. 

3.2 Singular Splitting Kernel. 

For the "twin" model ||], we have P(z) = 5{z - 1/2), P = I /2 s . We shall 
proceed here as before, by expanding the exponential in (|l2j), 



Us,t) = e-*± ^(^t) (47) 



n=0 

<fe(M) = 2(1- e ^)£i|L(_L) (48) 

n=0 

n-1 OO /nAfe 



A-' 

m=l fc=0 



1 



fc+m 



and applying the inverse Mellin transform term by term to find, 

= e-'f^(,-i) (50) 

n=0 ' ^ ' 
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d 2 (x,t) = 2(1 - e"') f] (s - -U (51) 

n=0 

n—1 oo (c).\k 

d n (x, t) = e("- 2 )*(l - e-^ 1 E ^" 1 2 m E ^TT X 

m=l fc=0 

^(^-2^) (52) 
For the factorial moment densities we have, 

KM) = .-E^*(.-sr)=^ £ ^ («) 

n=0 x ra=[— log 2 a;] 



C/ 2 (*,i) = 2(1 - e"*) E ^ ^ ' 



n=0 



OO 



= 2(l-e-*) E (54) 
n=[-log 2 a]-l 

u n ( x ,t) = ^-^(i-e-^-i^cr^E^rr^^ 



, n /c! V 2 fc + r 

m=l fc=0 
n-1 x ^ 

"¥" 



m=l fc=[— log 2 x]-m 

where [z] is the smallest integer n such that n > z and we denoted H the 
right continuous Heaviside step function. 
We shall use now the fact that, 

7 being an incomplete gamma function [26], so that, 
_ f 7 ([-log 2 x],2t) 

" 6 r([-io g2 *]) (5b) 

17 2 (M) = 2e (1-e ) r ([-log 2 x] - 1) (5?) 

U n (x,t) = e "*(l - e-)- 1 E Cfr 1 ^-^-™'^ (58) 

r([-log 2 x] -m) 
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3.2.1 Large t. 



For fixed x and t — ► oo the ratio of gamma functions goes to 1 and we have 
KNO scaling, U q /Uf — ► g!, by definition of the coefficients C^ -1 , as in the 
previous case. 



3.2.2 Small x. 

The functions t/q have a stairway shape superimposed to the general profile 
obtained by eliminating the integer part from the arguments. We shall work 
with this smoothed form in what follows. The density U q {Ui) is unaltered 
by these device, as can be explicitly checked, and it is in this density that 
we are interested. 

We shall now consider the form of these functions for large y = — log 2 (x) . 
To that end we use the following asymptotic expansion for large values of 
the parameter, obtained by a saddle-point approximation, 



j(y,x 



x y e 



1 + 



y-i 

r(y) = V^e-Vyy-^ 2 



1 x(x - 2) 
-+ o +••• 



y v 
l l 

+ 12y + 288y 2 + 



(59) 
(60) 



This way, we obtain, to first order in l/y, 



log 2 C/i = - ylog 2 (y) + y log 2 (2re) - - log 2 (y) + log 2 



+ 



24-1/2 



-t \ 



2vr/ 



+ 



(61) 



log 2 U 2 = -y log 2 {y) + y log 2 (2te) + - log 2 (y) + log 2 



2t - 1/2 ^ 

+ — + o 

y 



2TTt 



+ 



(62) 



The case of U n is slightly more complicated since it is a sum of several terms. 
Assuming y » 2t + n we have, for k < n, 



-f(y-k,2t) T(y-k+l) ^ y-k 
T(y-k) 7 (y-fc+l,2t) 2* 



> 1 



(63) 
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and then we can retain only the leading and next-to-leading terms n = n — 1 , 
n — 2, so that, 



U n * e nt (l - e 



-t\n-l 



■i 7(y~ w + l,2t) 
T(y-n+l) 



, (n-l){n-2) 2n _ 2 l (y-n + 2,2t) + 



T(y-n + 2) 



(64) 



We shall use a first order approximation for gamma functions in the first 
term and zeroth order for the ones in the second. Up to higher order terms 
we finally have, 



\og 2 U n = -ylog 2 y + ylog 2 (2ie) + (n- ^ log 2 y + est. + 

n(n — 3)~ 



1 

+ - 

y 



13 (n-1 n-2) 

2t + - -2t 

12 4 



(65) 



Proceeding as in the previous case, we arrive at the following expressions 
for the intermittency indices, 



1 + 



y (Hy) - M2t)) 

n- 1/2 



(4 



(n- 1) + 



y(ln(y) - In(2t)) 



+ 0(- 



(66) 
(67) 



Taking into account that y = — log 2 (x), we see that for — log 2 3> 2t + n we 
have the same spectrum of asymptotic intermittency indices as in the case of 
a uniform P(x), with logarithmic corrections independent of t to first order. 



4 Final Remarks. 

In the previous sections we considered two particular cases of one-species 
branching processes of the pure-birth type with binary fission, characterized 
by a uniform and a Dirac-delta splitting kernel, respectively. We exploited 
the fact that for these special kernels explicit solutions to the evolution equa- 
tions for inclusive distributions can be found. We have shown that for a fixed 
resolution in x, and t — > oo a KNO scaling regime is reached, while for t 
fixed and x — > factorial moments grow without bound, giving asymptotic 
intermittency indices ip q = q — 1 in both cases. Leading corrections to this 
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behavior for finite x are proportional to 1/lnx, and independent of t. Fur- 
thermore, they are positive, i.e., the asymptotic linear region is approached 
from above. 

The intermittent regime is defined as y2t -C \f—\nx for the uniform and 
as 2t <C — log 2 x for the singular kernel, for the first few indices with n ~ 1. 
It follows that x <C {x)(t) is a necessary condition, which seems difficult 
to fulfill in practical applications such as numerical simulations. In [28] a 
Monte Carlo experiment with three branching models is reported, of which 
only the "twin" model has been considered here. Their results, as shown in 
Figure 2 of pq] , are not inconsistent with our prediction for n < 4, within 
numerical uncertainties in their computation, although it should be noticed 
that for q > 4 their indices are systematically larger than q — 1, which is the 
maximum possible value [flTf ]. 

Even though our approach does not appear to be easy to extend to 
more complicated functional forms of the kernel P(x), we conjecture that 
the spectrum of intermittency indices is independent of it within the class 
of models considered in this paper. Our results may be useful to tune up 
numerical studies of these models, along the lines of [EH]. We believe that 
a richer intermittent behavior, as observed in experimental data, may be 
obtained with infrared divergent kernels, about which results have already 
been obtained [14]. 
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Appendix. 

We present here some basic equations and definitions from the theory of 
branching processes as used in sections 2 and 3. 

The starting point is the generating functional of the process, 4>[W, t; f] 
|20| , |2lf| , which is a function of the evolution parameter t and of the energy W 
of the initial particle, and a functional of the dummy function f{oj). Prob- 
ability densities will be obtained by functional differentiation with respect 
to /. Notice that u has dimension of energy, like W. 

The evolution of <f> is given by the equations [20, pif| , 



Ot 



[W,t;f] = J dzP{z)(^[zW,t-J]4>[{l-z)W,t-J}-(j){W,t-J^{kl) 

x [f{ziW)f{{\ - z)£W) - mw))j (A2) 

The first being the backward and the second the forward Kolmogorov evolu- 
tion equations [29]. These two equations have the same set of solutions [29], 
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and only the forward equation is used above. We use the initial condition 
(j)[W,t = 0; /] = f(W), corresponding to one particle in the initial state. 
Inclusive probability densities are defined as, 



D n ( Xl ,...,x n ;t) = W n 



5 n 0[W,t;f] 



5f(x n W)...5f( Xl W) 



(A3) 



/=i 



where we have already used the fact that since the splitting kernel P{x) 
depends only on energy fraction, and not on energy, the process is scale 
invariant and probability densities depend only on x\, . . . , x n [21|. Further- 
more, we see from this definition that D n is completely symmetric in its 
arguments. 

From the evolution equations for <f> we can derive the corresponding 
equations for D n . The forward one is given by (||) and the backward equation 
is, 



dDn 

dt 



\X\ , . . . , x n , i) 



\7r(n,K) 



(1-z) 



n—k 



D n —k 



.L : 



\ 



(A4) 



where the sum runs over all possible partitions vr(n, k) of (1, ... , n) into two 
sets . . . ,ik), (ik+i, ■ ■ ■ ,in) such that i\ > %2 > • • • > ik and i^+i > • • • > 
i n , with k = 0, . . . , n; and we used the relation, 



W 



[zW,t;f] 



5f(xW) 



1 -D X ( X -,t 

z \z 



(A5) 



and its generalizations to higher order densities D n . 
Analogously, exclusive distributions are defined by, 



E n ( Xl ,...,x n ;t) = W n 



5 n <t>[W,t-f] 



5f(x n W)...5f(x 1 W) 



(A6) 



/=o 



and satisfy, 
8E„ 



dt 



(xi , . . . , x n , i) 
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E k ( — , . . . , — ;t] E n _ k ( — — , . . . , — ;t ) — E n (x\, ... ,x n ;t) 



z z 



\ 

;a7) 
I 



^r(x\, . . . ,x n ,t) = -nE n (xi, . . . ,x n ,t) + - V P\ — — — )x 
dt n k ^ =1 \x J+ x k J 

x . Jj _ x E n-i(xi, ■ ■ ■ ,Xj + x k , . . . ,xj,, . . . ,x n ,t) (A8) 

j 

The initial condition is given by E n (x\, . . . ,x n ; t = 0) = 5i n 5{xi — 1). 
The relation between D n and E n may be explicitly written as, 



D n ( Xl ,...,x n ;t) = N(N-l)---(N-n + l) 



oo 

N(N -1)---(N -n + 1) x 

N=n 
»1 

dy N= i ■ ■ ■ dy N E n (xi , . . . ,x n , yx, . . . , y n ; t) (A9) 

o 

where, once again, we have made use of the symmetry with respect to permu- 
tations of the arguments of these distributions. Notice that our definition 
of E n is slightly different from that of ^l), since we normalize E n to the 
multiplicity distribution, 

/ dxi ■ ■ -dx n E n (xi, . . . ,x n ;t) = P n (t) (A10) 
Jo 

from whence the normalization of D n to factorial moments results. 

It is not difficult to convince oneself that solutions to the evolution equa- 
tions have the form, 

E n (x 1} ...,x n ;t) = P n (t)K n (xi, ...,x n ) (All) 

where K n is defined by this equation. The evolution equations for P n (t), the 
multiplicity distribution, can be obtained by integration of the equations for 
E n . Since these are well-known [11, 23 1, we shall quote only the solutions, 

P„ (t) = i ( w^ir 1 (Ai2) 



(n) V (n) 

(n) = e* (A13) 

From this expression, and the evolution equations for E n , recurrence rela- 
tions (backward and forward) for K n (xi, . . . , x n ) can be found. We shall 
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only write the forward one, 

2 n ( x \ 1 
K n (xi,...,x n ) = — 22 M — T — — T x 

n ( n - 1 ) k ^~L 1 \Xj+XkJXj+X k 

K n _i(xi, ...,Xj + x k ,...,Xk,...,x n ) (A14) 

3 

with K\(x) = 5{x — 1). From this relation and P(x) = P(l — x) it follows 
that, 

f 1 1 

/ dx 1 ---dx n x 1 K n (xi,...,x n ) = - (A15) 

Jo n 

and then, 

00 ^ ln(n) 



(x) = V / dxj • • • cfe„ xi£' n (xi, . . . ,x n ;t) 

n=l 



(n)-l 



(A16) 



independently of the detailed form of P[x). 

In fl2]| , it is shown that a branching process satisfying, 

n I n \ 

K(xi,...,^„;t) = 11/(3*;*)* h-X^y ( A17 ) 
i=i v i=i / 

is not intermittent. From the factorization property discussed above, we see 
that such a branching process does not exist, at least within the class of 
models considered here. 
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